Precessing Vortex Motion and Instability in a Rotating Column of Superfluid 3 He-B 
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The flow of quantized vortex lines in superfluid 3 He-B is laminar at high temperatures, but below 
0.6 T c turbulence becomes possible, owing to the rapidly decreasing mutual friction damping. In 
the turbulent regime a vortex evolving in applied flow may become unstable, create new vortices, 
and start turbulence. We monitor this single- vortex instability with NMR techniques in a rotating 
cylinder. Close to the onset temperature of turbulence, an oscillating component in NMR absorption 
has been observed, while the instability generates new vortices at a low rate ~ lvortex/s, before 
turbulence sets in. By comparison to numerical calculations, we associate the oscillations with spiral 
vortex motion, when evolving vortices expand to rectilinear lines. 

PACS numbers: 67.57.Fg, 47.32.-y, 67.40.Vs 



In superfluid dynamics a longstanding goal has been 
to account for the formation of all vortices. A quantized 
vortex line is a topologically stable structure of a coher- 
ent order parameter field, with fixed circulation of su- 
perflow. What is then the mechanism by which sudden 
burst-like turbulent vortex proliferation is started from 
one single vortex loop which is evolving in applied flow? 
The common view holds that turbulence, once started, is 
sustained via loop formation and reconnections, as seen 
in well-developed thermal counterflow turbulence of su- 
perfluid 4 He-II pp. However, as emphasized by Schwarz 
[2], to start turbulence, for instance in linear flow along 
a circular pipe, vortices have to be pinned in "vortex- 
mill" configurations which continuously inject new vor- 
tex loops downstream in the applied flow. Nevertheless, a 
turbulent burst of vortex formation is observed to evolve 
from a single seed vortex loop in rotating 3 He-B, with no 
permanently pinned vortices [3]. 

A recent explanation [4] of this controversy concludes 
that a vortex evolving in applied flow may become unsta- 
ble while interacting and reconnecting with the container 
wall. In 3 He-B this single-vortex instability becomes pos- 
sible with decreasing temperature below an onset temper- 
ature T on , where mutual friction damping has dropped to 
sufficiently low value: a(T) < 1. In the rotating cylinder 
the instability leads to a sudden transition to the equi- 
librium vortex state. This process is most conveniently 
studied in the onset temperature regime T ~ T on , by 
monitoring the number of vortex lines as a function of 
time at constant rotation velocity, as a response to the 
injection of a seed vortex [5]. The events following the in- 
stability have been reviewed in Refs. [7] for the case of 
a long rotating column. Here the end point of an evolv- 
ing seed vortex describes a spiral trajectory along the 
cylinder wall while it expands towards its stable state as 
rectilinear vortex line. In the onset region T ~ T on , the 
spiral vortex motion is sometimes observed to give rise to 
an oscillating NMR absorption signal which is examined 
in this report. 




FIG. 1: Numerical calculation of vortex evolution in a ro- 
tating cylinder of 3 He-B: Rotation is suddenly increased at 
t = from ft « 0.03rad/s to Q { = 0.2rad/s. There are 22 
vortices in this sample, of which two in the outermost ring 
(lying opposite to each other) have been initially bent to the 
cylindrical wall, to break cylindrical symmetry. In the later 
snapshots at Qf, the two short vortices expand towards the 
top and bottom end plates of the cylinder, while other vor- 
tices have contracted to a central vortex cluster. Parameters: 
R 3 mm, L 30 mm, P 29.0 bar, and T 0.4 T c (which 
corresponds to a = 0.18 and a = 0.16 [5]). 

Numerical illustrations: To visualize the motion of 
evolving vortices in rotating flow, a calculation is pre- 
sented in Fig. [I] using the vortex filament method de- 
scribed in Ref. [5]. The initial configuration in Fig. [I] 
with two curved vortices which connect at one end to the 
cylindrical side wall, mimics the equilibrium vortex state 
in a real rotating experiment. Depending on the rotation 
velocity Q[ and the residual misalignment between the 
rotation and the sample cylinder axes (< 1° in the setup 
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FIG. 2: Comparison of azimuthal (top) and longitudinal 
vl z (bottom) velocities when different numbers of 2iVb vor- 
tices expand in the setup of Fig. [I] With 22 vortices in total, 
1, 2, 3, or 9 vortices of equal initial length 0.4 L have been bent 
both on the top and bottom to the cylindrical wall, to mimic a 
tilted cylinder, as shown in Fig.[l] The spiralling motion both 
upwards and downwards along the rotating column is calcu- 
lated as a function of time. The average velocity of the vor- 
tex ends on the cylindrical wall in the upward and downward 
moving bundles is plotted during the time needed to reach 
the respective end plate. The actual velocity is the number 
on the vertical scale times av(Qf , R, N) or (1 — a')v(Q{ , i?, iV), 
where N = 22 - N h . The case iV b = 22 is different: here all 
22 vortices are initially bent at one cylinder end to the wall 
and, after the rotation increase to Qf, a vortex front starts 
expanding towards the other vortex-free end of the cylinder. 
This situation is examined further in Fig. [3] 

of Fig. |5j, a certain fraction of the peripheral vortices 
ends on the cylindrical wall [9] . 

In Fig. [I] rotation is increased in step-like manner at 
t = from the equilibrium vortex state at (]j / to a 
final stable value Qf. The N = 20 rectilinear vortices 
are thereby compressed to a central cluster with an areal 
density n v = 2Qf /k by the surrounding counterflow (cf) 
at the velocity v = v n — v S: the difference between the 
velocities of the normal and superfluid fractions. The 
highly viscous normal component we consider to be in 
solid-body rotation: v n = Qxr. Outside the compressed 
cluster the cf velocity is then given by 



v(Qf, r, N) = v n — v s = Qfr 



2ttt ' 



(i) 



where ft = h/(2ms) is the quantum of circulation. The 
two short vortices now expand in this cf in spiral motion 
towards the top and bottom end plates, respectively. The 
velocity vl of a vortex line element with tangent s is 
obtained from the equation of motion [10] 

vl = v s + as x (v n - v s ) - as x [s x (v n - v s )] , (2) 



where a is the dissipative and a' the reactive mutual 
friction coefficient. Since the vortex end is perpendicular 
to the cylindrical wall, it has from Eq. ^ a longitudinal 
velocity v^ z = av(Cl{, R, N) and an azimuthal component 
= (1 — a')v(Qf , i?, N). In the rotating coordinate 
system, points in the opposite direction from the 
rotation of the container. The resulting end point motion 
gives us a simple picture of the spiral trajectories along 
the cylindrical wall, although other parts of the vortex 
also contribute to its motion, in particular its curvature 
where it connects to the cylindrical wall. During their 
expansion the two evolving vortices are wound around 
the central vortex cluster as a helix whose pitch depends 
on the cf velocity i>(Qf , R, N). Ultimately, the helix will 
unwind, when the vortex ends slip along the flat end 
plates of the cylinder, so that the final state is composed 
of only rectilinear vortex lines. However, since the cf 
velocity is close to zero at the border of the cluster, the 
unwinding is a slow process. 

The calculated velocities of the vortex ends in Fig.[l]are 
v Lz w 0.84 aQR w 0.96 av(Vt u R, N) and v L(f) w 0.73(1 - 
a')QR w 0.83(1 - a') v(Q f ,R,N). The wave length of 
the spiral trajectory is thus A = 2i\Rv\, z lv\,^ ~ 5 mm 
and the period p — 2itR/v\ J( f ) « 50 s. In Fig. [2] the spiral 
vortex motion is analyzed further for the same experi- 
mental setup as in Fig. [T] but when several vortices are 
expanding simultaneously. The initial starting state for 
these calculations is one where a specified number (iVb) 
of nearest-neighbor vortices is bent to the cylindrical wall 
at the same distance of 0.4 L from the end plate both at 
the bottom and top of the cylinder, similar to the config- 
uration in Fig. [l] (at t = 0). When rotation is increased 
from Q[ to fif, one set of curved vortices starts ex- 
panding towards the top and one towards the bottom end 
plate. The axial (bottom panel) and azimuthal velocities 
(top panel) have been plotted as a function of time, for 
different numbers of 2N\> vortices in spiral expansion. 

Two observations can be made from Fig. [2j The vor- 
tices, which expand in one direction and initially start off 
in close proximity of each other, remain a close bundle 
during their spiral motion. Vortex bundling is a general 
phenomenon which gives rise to the formation of larger 
eddies, similar to the eddies in classical viscous fluid mo- 
tion. Here the bundles do not disperse, but are well pre- 
served especially during the later part of the expansion 
in a state of steady propagation, since all N\> vortices 
travel roughly at the same speed both axially and az- 
imuthally along the cylinder. Thus the spread in veloc- 
ities among the different vortices within the bundle and 
also the changes in the velocities as a function of time 
remain relatively small. The second note about Fig. [2] is 
that both vl z and are rather insensitive to the num- 
ber of vortices spiralling in the bundle. This is because 
the total number of vortices N = 22 corresponds here 
to a very low equilibrium rotation velocity ~ 0.03rad/s 
compared to the actual rotation velocity of 0.2 rad/s. For 
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FIG. 3: Calculated azimuthal front velocity Vf^ (in the rotat- 
ing coordinate system) , as a function of the number of vortices 
Nf which compose the front. No prior vortex cluster exists 
in the center of the cylinder. The zero temperature limit of 
Eq. [3] (solid curve) is compared to calculations at finite tem- 
peratures (data points). The calculations with Nf < N eq are 
for R — 3 mm, L — 30 mm, il — 0.2 rad/s, where iV eq ~ 170 
vortices: (o) a = 0.1, a = 0; (x) a = 0.18, a = 0.16, 
T w 0.4 T c . The two calculations with Nf ~ iV eq are for 
R — 1.5 mm, L — 40 mm, Q — 1 rad/s, where N eq ~ 210 vor- 
tices: (y) ot = 0.040, ot = 0.030, T « 0.3 T c ; (A) a = 0.18, 
a — 0.16, T w 0.4 T c . On the vertical scale the normalized 
velocity Vf^/[(1 — a')QR] is plotted. Vf^ is taken as the front 
velocity when a stable value has been reached, which gener- 
ally happens when the front has propagated 2/3 of the length 
of the cylinder. 

better comparison with measurements the final rotation 
velocity fif should be increased even more, which calls 
for time-consuming calculations. 

A special case is that where all existing vortices are re- 
leased simultaneously from one end of the rotating cylin- 
der and expand as a vortex front towards the vortex-free 
end of the cylinder {i.e. where no prior central vortex 
cluster exists). In Fig. [2] this case is the example with 
iVb = 22. Such a calculation models the experimental 
situation after a turbulent burst, like in Fig. [6] (although 
there the burst occurs higher up in the column, which 
then starts both an upward and downward moving front). 
At temperatures below about 0.45 T c the front, once it is 
formed from Nf vortices, maintains its constant narrow 
width during the spiralling propagation along the col- 
umn [6j [7] . The longitudinal propagation velocity Vf z of 
the front has been measured and analyzed in Ref. [TT] . 
The azimuthal velocity Vf<j> (in laboratory coordinates) 
is obtained in the continuum model at zero temperature 
(where friction vanishes) from the equation 

^7V f l/2 1n(7V eq /7V f ) + l/4 

f0 2KR 1 - Nf /(2JVeq) ' [ > 

where the approximate continuum value for the number 
of vortices in the equilibrium state is N eq = ttR 2 (2CI/ 'k). 
When Nf — > N eq , this expression gives Vf^ = \ Q.R. This 
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FIG. 4: (Left) Rotating vortex cluster from Fig.[l]at t — 60s, 
viewed at an angle deviating by 3° from vertical. Only the 
central part of the sample cross section in the x - y plane is 
shown. Helical Kelvin waves are induced on the vortex lines 
in the cluster by the two vortices spiralling outside the cluster. 
In the early stage of the expansion, the waves start at roughly 
half- height of the cylinder, i.e. from the points where the two 
short vortices merge with the cluster (see Fig. [I). (Right) Side 
view through the vortex cluster (t = 60 s) . A vortex spiralling 
around the cluster induces radial displacements on the out- 
ermost vortices within the cluster. The displacements start 
longitudinally propagating Kelvin waves on these vortices. 

limit was discussed in Ref. [12] . In Fig. [3] we compare 
Eq. |3| to numerical calculations at finite temperatures 
and finite friction for a few examples. The results for 
the normalized azimuthal front velocity Vf</>/[(! — a')QR] 
drop below the estimate from Eq. ([3]), but by less than 
10%. Thus also at finite friction Vf^ « \{1 — a f )QR is a 
good approximation, which applies when the number of 
vortices in the front is large (TV < A^ eq ). 

While spiralling around the cluster in Fig. [TJ the two 
curved vortices may reconnect with each other or with 
the outermost lines in the cluster. These reconnections 
do not increase the vortex number. As can be seen from 
Figs, [l] and |4j a spiralling vortex distorts the cluster by 
inducing propagating helical Kelvin waves on the vortex 
lines, which results in oscillations of the cluster around 
its equilibrium position [13] . The amplitude of the waves 
is comparable to the inter- vortex distance in the cluster, 
~ 1/ ^Jn^ 1 and initially the wavelength and period of the 
excited Kelvin waves have similar values as the precessing 
vortex motion outside the cluster. 

At high friction (a > 0.18) the number of vortices re- 
mains strictly constant in Figs. [T] — [4] At low friction 
(a < 0.1) the calculations on the precessing vortex front 
in Fig. [3] may display some increase in vortex number, 
when N < 7V eq . Figs. [T]- [2] thus illustrate the situation 
above T on where no increase in vortex number occurs. 
Experimentally T on depends on the applied flow velocity 
(~ [Qf — Q[]R) and on the number (Nf) and initial con- 
figuration of evolving vortices [5j [3 [14]. In most cases 
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T on is found to be higher than the temperature of the 
present calculations, 0.40 T c . The reason for this differ- 
ence appears to be the high stability of rotating flow in 
a circular cylinder, which according to our calculations 
[5] appears to be a particular characteristic of a cylinder 
with a circular cross section and ideal wall properties, 
i.e. without surface friction or pinning. The most likely 
section of an evolving vortex to become unstable is the 
curved piece ending on the cylindrical wall in the max- 
imum flow u(f2f, i?, TV), while the most stable parts are 
the sections which reside in low cf close to the central 
vortex cluster. The extreme case are rectilinear vortex 
lines in the cluster which are experimentally found to re- 
main stable even in sinusoidally modulated rotation at 
temperatures > 0.3 T c . 

Experiment: The sample is cooled in a rotating nu- 
clear demagnetization cryostat at 29 bar liquid pressure 
within a fuzed quartz cylinder, with radius R = 3 mm 
and length L = 110 mm (Fig.|5|. The number of vortices 
is measured non- invasive ly simultaneously at both ends 
of the long cylinder, by monitoring the NMR absorption 
line shape. Both NMR detectors consist of circular split- 
pair coils, with a separation of 9 mm between the coil 
halves and with their common axis aligned transverse to 
the cylinder. The coils are wound from superconduct- 
ing wire and are part of a LC resonator with a Q value 
of - 1 • 10 4 , coupled to a GaAs MESFET preamplifier 
located within the cryostat at 4.2 K temperature. Both 
NMR spectrometers operate independently at constant 
frequency in an axially oriented polarization magnetic 
field with a linear sweep. The NMR setup is described 
in more detail in Refs. [3l [T5] 

With vortex-free cf, the NMR absorption line shape in- 
cludes a prominent cf peak which is shifted far from the 
Larmor frequency, while in the equilibrium vortex state 



the absorption drops to zero at the location of the cf peak 
[T6] . The height of the cf peak increases with the mag- 
nitude of the azimuthal cf velocity v(Q,R,N) and thus 
decreases with the number of vortices N in the central 
cluster. The correspondence between peak height and 
vortex number N can be obtained from calculations of 
the order parameter texture or from measurements. The 
measurements which we discuss below use remanent vor- 
tices [14J as the seeds which evolve in the applied flow 
after the rotation increase to fif . Such a measurement is 
performed at constant temperature, by first decelerating 
rotation to zero from a state with a large number of vor- 
tices. Zero rotation is then maintained for a period At, to 
allow vortices to annihilate, with the exception of a few 
remaining dynamic remnants [H] . The final step is to 
increase Q at a fixed rate (typically dd/dt < 0.02rad/s 2 ) 
to fif, where it is kept constant while the buildup in 
the number of vortex lines N(t) in the central cluster is 
recorded as a function of time. 

Single vortex instability: Well above the onset tem- 
perature, the cf peak height settles at f2f to a stable value 
and remains constant. Well below onset the cf peak 
collapses to zero when the equilibrium vortex state is 
formed. Around the onset temperature T « T on , both 
types of behavior occur randomly. If the equilibrium 
vortex state is here formed, then in perhaps one third 
of the cases the collapse of the cf peak can be preceded 
by an initial slow decrease in peak height. The slow re- 
duction in peak height corresponds to slow generation of 
new vortices. It is the signature of the single- vortex in- 
stability, as described and analyzed in Ref. [5]. It can 
only be observed in the onset regime T ~ T on ; at lower 
temperatures the instability proceeds too rapidly to be 
monitored with our techniques. The slow vortex genera- 
tion is terminated in a turbulent burst which takes place 
as a localized event in some short section (of length ~ R) 
of the long sample column. In this burst enough vortices 
are created to start the evolution towards the equilibrium 
vortex state. From the site of the burst the vortices prop- 
agate in spiral motion as a front both up and down along 
the cylinder, leaving behind a vortex bundle composed 
of helically twisted vortices [6j [7] . 

The NMR absorption responses of rapid and slow tran- 
sitions to the equilibrium vortex state are compared in 
Figs. [6] and [7| respectively. These were recorded in two 
consecutive measuring runs at the same temperature us- 
ing the same measuring procedure. The order of the 
two measurements was such that Fig. [7] was measured 
first and Fig. [6] next. The responses prove to be differ- 
ent owing to the different value of the rotation veloc- 
ity Qf, which in turn is caused by the stochastic na- 
ture of the remanent vortex injection process: in the 
former case the waiting period at zero rotation prior to 
the measurements was At = 40min, while in the lat- 
ter case it was 22 min. Thus in Fig. [6] the remanent 
vortices are expected to be fewer and smaller loops. In 
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FIG. 6: NMR response during a rapid transition from vortex- 
free flow to the equilibrium vortex state at constant externally 
controlled conditions. The top and bottom detectors monitor 
the NMR absorption at different values of constant NMR field. 
The cf peak height is recorded with the top detector, while 
the bottom detector monitors the peak close to the Larmor 
edge. Time t — marks the moment when Qf — 1.61rad/s 
is reached. The time delay between the sudden rise in the 
bottom response (at t = 40 s) and the start of the rapid col- 
lapse in the top (at t = 78 s) is caused by the longitudinal 
motion of the two vortex fronts along the column, controlled 
by a = 0.33 0. 

both cases the sample was then accelerated to rotation at 
Cl = 0.004 rad/s 2 . The acceleration was interrupted oc- 
casionally, to check at constant rotation whether vortex- 
free rotation continued to persist over long periods in 
time. In Fig. [6] (Fig. [7]) rotation could be increased until 
1.61 rad/s (0.90 rad/s), before any indication of vortex 
formation was noticed. Time t = coincides with the 
moment when Qf = 1.61 rad/s (0.90rad/s) was reached. 
After 40 s (170 s) the first response from vortices is ob- 
served in the bottom (top) spectrometer. This should 
be compared with the time it takes for a small vortex 
loop to expand in vortex-free flow at fif from one end 
of the cylinder to the other, to become a rectilinear line, 
« L/(aSlfR) = 70s (120s). In Fig. [7] the start of vortex 
formation is thus unusual because of the long 170 s delay, 
before the generation of new vortices starts. The expla- 
nation for the long delay must be associated with the 
properties of a small remanent vortex loop with a radius 
close to the critical value of ~ 4 /im needed to overcome 
the barrier for spontaneous expansion [T7] . 

Counterflow peak response : In Figs. [6] and [7] the 

peak heights of two different absorption maxima of the 
NMR spectrum are recorded as a function of time at 
constant externally controlled conditions. These are: (1) 
The cf peak height, which initially in vortex- free flow is at 
maximum and zero in the final equilibrium vortex state. 
(2) The Larmor peak height, which initially is close to 
zero, but in the final state different from zero. While the 
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FIG. 7: Precursory vortex formation which terminates in a 
sudden turbulent burst. This is a repetition of the mea- 
surement in Fig. [6] but at a lower rotation velocity Qf = 
0.90 rad/s. The difference is the prolonged slow vortex gener- 
ation via the single- vortex instability, denoted by the dashed 
line. Superimposed comes the large- amplitude quasi-coherent 
oscillation in the cf peak height at a frequency which is re- 
lated to the azimuthal vortex motion. The latter is controlled 
by I -a' = 0.76 0. 

Larmor peak response is qualitatively similar in the two 
examples, the drop in the cf peak height happens dif- 
ferently. Fig. [6] shows the generic response from a rapid 
transition, while in Fig. [7] the decay in the cf peak height 
starts first with a more gentle linear decrease (dashed 
line) where the number of vortices gradually increases 
within the top detector, as generated by the single- vortex 
instability. This continues for more than 100 s, before the 
rapid peak height decay starts. (Also in Fig. [6] at higher 
Qf a short vestige of slower peak height reduction can 
be distinguished, but the rate is faster and it only lasts 
for ~ 4 s.) However most importantly, in Fig. [7] the slow 
decrease in peak height differs from usual examples in 
that the cf peak height also displays well resolved oscil- 
lations of large amplitude. The characterization of these 
oscillations is the central issue in this report. 

According to the measured calibration, the linear de- 
crease in the cf peak height, as marked by the dashed 
line in Fig. represents a vortex formation rate N w 
1.2 vortices/s. The later rapid collapse is caused by the 
arrival of the vortex front [TT] to the bottom end of the 
top coil and its subsequent travel through the coil. At the 
arrival of the front, the central cluster within the top coil 
contains ~ 140 vortices (to be compared with N eq ~ 840 
in the equilibrium vortex state). The front incorporates 
approximately all the remaining vortices needed to fill 
the sample with the equilibrium vortex state. The front, 
composed of these additional vortices, travels in spiral 
motion along the cylinder around the central cluster. Be- 
hind the front the expanding vortices are wound in a heli- 
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cally twisted configuration which later slowly relaxes [12] . 
At that point the macroscopic cf, which was generated 
by the rotation increase to Qf, is completely removed. 

Larmor peak response: The second signal trace in 
Figs. [6] and [7] is recorded by scanning the absorption max- 
imum close to the Larmor edge with a rapidly moving 
field sweep. Here the sudden sharp absorption increase 
signals the arrival of the vortex front, as it passes through 
the top edge of the bottom detector. The increasing ab- 
sorption is generated by the axially flowing supercurrent 
which is produced by the helically twisted vortices behind 
the front [18 . The subsequent exponential decay of this 
absorption is caused by the unwinding of the twist while 
the vortex ends slip along the top and bottom end plates 
of the sample cylinder and the vortices are converted to 
rectilinear lines [T2] . Thus the Larmor absorption re- 
sponse forms a transient peak, followed by a stable final 
absorption level, which is representative of the equilib- 
rium vortex state, after the twist has relaxed [6]. 

In Fig. [7] a delay time t& ~ 144 s separates the mo- 
ments when the vortex fronts moving up and down along 
the column pass through the bottom edge of the top 
detector and the top edge of the bottom detector, re- 
spectively. This distance is 82 mm. The longitudinal 
propagation velocity Vf z of the front has been measured 
to be Vi z / (SIR) ~ 0.28 at 0.45 T c [IT], in a situation 
when there is no prior central vortex cluster. Assuming 
that the propagation velocity of the front in the presence 
of a central vortex cluster is ~ v(Clf,R,N) V{ Z /(CIR), 
we estimate that, with a vortex formation rate of N ~ 
1.2 vortices/s, on an average TV « 230 vortices populate 
the central cluster at the time when the front travels from 
the location of the turbulent burst to the top edge of the 
bottom detector. The delay time of 144 s can then be 
used to localize the site and the moment of the turbulent 
burst. Here the burst proves to happen a couple of mm 
below the top detector at t « 290 s, or some five seconds 
before the absorption response in the top detector starts 
to collapse. In Fig. [6] at a higher flow velocity all events 
occur faster: with a delay of 33 s between the first bottom 
and top responses, the burst happens at t = 28 s some 
19 mm above the top edge of the bottom coil. 

Spiral vortex motion: The generic signal from slow 
vortex generation in the cf peak height is a relatively 
smooth decay, while the many periods of well-resolved 
quasi-coherent oscillation in Fig. [7] are unusual. If both 
the top and bottom detectors are tuned to monitor the 
cf peak height simultaneously, then the two recordings of 
the oscillations look similar, but neither the amplitudes 
or the frequencies are exactly identical. Since the distri- 
bution and number of evolving vortices varies along the 
column, this is expected. Obviously the oscillating signal 
is driven by the spiral vortex motion, which is the only 
source for precession at a frequency around 0.1 Hz. In 
Fig. [8] we have extracted the difference in time between 
the maxima and minima in the oscillation of Fig. [7| This 
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FIG. 8: Time difference between the maxima and minima in 
the oscillating cf peak height in Fig. [7| The line is obtained 
by inserting the appropriate values from Fig. [7] in Eq. Q. 

data set consists of 13 full periods. The average period 
is 9.2 s, but the period also appears to be increasing with 
time. The spiral vortex motion is expected to have the 
frequency 



= e (1 - a')v(Slf , R, N)/(2ttR) , 



(4) 



so that oc —N. We assume that the azimuthal vor- 
tex velocity is of the form = e(l — a')v(Qf, R, iV), 
where s < 1 is a numerical factor, which in the case of 
Fig. [2] is about 0.8. The measuring setup in Fig. [5] incor- 
porates also an instrument factor: The sensitivity of the 
detector coil pair is not constant across the cross section 
of the sample, but increases towards the windings. The 
cf peak height is regulated by the azimuthal vortex-free 
flow, i.e. the height depends on the velocity v(Qf,R,N) 
at the cylindrical wall. This doubles the frequency of 
an azimuthally precessing asymmetry in the cf velocity 
v(Qf,R,N). For simplicity we assume that e = \ which 
then compensates for the frequency doubling owing to 
the sensitivity pattern. 

If we take N ~ 1.2 vortices/s and N — at the moment 
when vortex formation starts (at t = 170 s) in Fig. [7| we 
obtain the line through the data in Fig. [8] The best 
fit would have a twice larger slope, but in view of the 
large scatter in the signal, which is only quasi-coherent, 
the agreement with the calculated line appears accept- 
able. Using this approach, we have extracted 1 — a' from 
seven runs with slow oscillatory cf peak decay. In these 
seven cases the single- vortex instability proceeds at a rate 
N ~ 0.4 — 2 vortices/s, while fif is in the range 0.8 — 
1.2rad/s. As seen in Fig. [9j the results agree with the 
measurements on 1 — a' (T) in Ref . [8] . 

Discussion: The oscillatory response in the cf peak 
height in Fig. [7] is a rare event where several precondi- 
tions seem to be fulfilled: 

1) New vortices have to be created at a slow rate 
N ~ 1 5 which requires that T w T on . 
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FIG. 9: Consistency test on the analysis of the cf peak height 
oscillations from seven different measuring runs. The reactive 
mutual friction coefficient 1 — a' has been extracted using 
Eq. Q and is compared to the measurements in Ref . [8] . The 
line is a fit from Ref. [8]. 

2) In all examples the longitudinal vortical transit time 
through the 9 mm long detector coil pair is roughly 10 s 
and equal to the oscillation period in the cf peak height. 

3) The oscillations start immediately when the slow 
vortex generation switches on. In Figs. [6] and [7] the reso- 
lution in the measurement of N is roughly one vortex, i. e. 
the reduction in cf peak height per one vortex is approx- 
imately equal to the amplifier noise (when N <C N eq ). 
This signal/noise resolution applies to the reduction in 
cf peak height, when one new vortex is placed in the cen- 
ter of the cylinder. Presumably the signal from a single 
vortex spiralling along the cylinder wall is larger. 

4) To explain the quasi-coherent oscillation, the spiral 
motion of several vortices has to be reasonably coherent: 
the 13 oscillations in Fig. [7] cannot be attributed to a 
single vortex since the transit time for a vortex end to 
pass through the detector coil pair is equal to one period. 
Three features may help to create and preserve coherence 
in spiral motion: (i) When a vortex end on the cylindrical 
wall creates a new loop, one end of the new loop starts 
spiralling close to the original end. (ii) Both detectors 
are close to one of the end plates of the cylinder and, as- 
suming that the single- vortex instability occurs randomly 
everywhere on the cylindrical wall [5 , almost all vortices 
approach the detector from the direction of the far end of 
the cylinder, (iii) In numerical calculations the evolving 
vortices, which start to spiral as a close bundle, tend to 
remain in a bunch while propagating along the cylinder. 

A complete explanation has not been worked out for 
the mechanisms which gives rise to the oscillating sig- 
nal in the cf peak height. To understand the coher- 
ence over many periods of spiral motion we need bet- 
ter agreement between vortex-dynamics calculations and 
the experiment, in particular concerning the generation 
of new vortices by the single-vortex instability in a ro- 



tating circular cylinder [5]. As to the amplitude of the 
oscillations, spiral vortex motion gives rise to oscillatory 
displacements of the order parameter texture from cylin- 
drical symmetry, as seen in Fig. gj To expl ain the signal 
from such oscillations, numerical calculations are needed 
on distortions of flare-out textures from cylindrical sym- 
metry and on the resulting NMR line shapes. 

Conclusions: We have studied the single- vortex in- 
stability of evolving vortices in the turbulent tempera- 
ture regime of 3 He-B. In applied flow the single- vortex 
instability generates new vortices and becomes the pre- 
cursor mechanism which starts turbulence. In a rotating 
circular cylinder the instability occurs while remanent 
vortices, for instance, expand in spiral motion towards 
their stable state as rectilinear vortex lines. The spiral 
motion has been examined here with numerical calcula- 
tions, but we also argued that oscillations in the NMR 
response bear direct evidence for the precessing motion. 
Other kinds of oscillating responses from precessing vor- 
tex motion have been reported before. A remarkable ex- 
ample is the unwinding of trapped circulation around a 
thin wire suspended along the axis of a cylindrical con- 
tainer, known as the Vinen vibrating wire experiment 
[T9] . The spiral motion of an evolving vortex around a 
central cluster of vortex lines provides another example 
of precessing vortex signals. 
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